Ion-mediated RNA structural collapse: effect of spatial confinement 
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Abstract 



PvNAs are negatively charged molecules residing in macromolecular crowding cellular environments. 
Macromolecular confinement can influence the ion effects in RNA folding. In this work, using the recently 
' go developed tightly bound ion model for ion fluctuation and correlation, we investigate the confinement effect 
on the ion-mediated RNA structural collapse for a simple model system. We found that, for both Na + and 
, O-f , Mg 2+ , ion efficiencies in mediating structural collapse/folding are significantly enhanced by the structural 
confinement. Such an enhancement in the ion efficiency is attributed to the decreased electrostatic free 
energy difference between the compact conformation ensemble and the (restricted) extended conformation 
ensemble due to the spatial restriction. 
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Introduction 



Cellular functions of nucleic acids are intrinsically related to their folding in cellular environments ([]]- 
3|). Since DNAs/RNAs are highly negatively charged molecules, folding into compact native structures 



requires metal ions to overcome the st rong Co ulombic repulsions (ll|)- (|13r) . Furthermore, RNAs in vivo are 
surrounded by many other molecules (|14l - ll6|) . and the volume percentage of macromolecules in cells can 
reach 40% (|14l - ll6[) . The presence of other macromolecules around an RNA can cause spatial confinement 
for DNA/RNAs folding. The macromolecular crowding may greatly influence the RNA folding (|l7l-l2ll). 
However, the problem of how the spatial constraint affects the ion-induced RNA folding is not fully solved. 
In this work, we will construct a simplified structural model to investigate the macromolecular crowding 
effect on ion-mediated RNA collapse/folding. 



For RNA secondary structure folding (tJIsl), the thermodynamic parameters at standard salt condition 
(namely, 1M Na + ), have been experimentally determined. Parameters for other ionic conditions have also 
been fitted from experimental data and theoretical calculations (|22|) - (|33l) . These parameters have led to 
many accurate predictions for RNA/DNA secondary structures and folding kinetics (|22|) - (|33l) . For tertiary 
structure folding, however, our understanding of the role of ions, especially the effects of multivalent ions 
such as Mg 2+ , remains incomplete (fl])- (fl3"|) . Experimental and computational studies showed that metal 
ions with higher charge density are more efficient in promoting RNA collapse and stronger ion-correlations 
could potentially contribute to the efficient role of Mg 2+ ions (|34t) -(l40l). Another series of experiments for 
short DNA helices suggested the existence of helix-helix attractive force at high multivalent ion concen- 
tration and repulsive force at low ion concentration 00. Furthermore, experiments on a paradigm 
system consisting of two helices tethered by a loop suggested that Mg 2+ is > 10-fold more efficient than the 
predictions from the Poisson-Boltzmann (PB) theory in promoting structural collapse (|43l 445|) . 

In the cells, the presence of other molecules can cause an crowded environment that can cause spatial 
confinement for the conformations of RNAs and consequently influence the free energy landscape of RNAs. 
For proteins, macromolecular crowding has been shown to play an important role in folding stability, struc- 
ture, and kinetics I46l-l51|). For RNAs, however, the study of macromolecular crowding effect on 



folding is relatively new (|52|) . especially on its impact on the role of ions in folding (|19h . 

RNA tertiary structure collapse involves high charge density, quantitative description for the role of 
metal ions is challenging, especially for multivalent ions. Molecular dynamics studies for ion-RNA in- 
teractions have provided novel insights into ion and structure/sequence- specific forces in RNA folding, 
especially the force arising from site-specific binding of ions (|53|) -(p8r). The two classic polyelectrolyte 
theories: the counterion condensation (CC) theory (59) and the PB theory (|60|) - (|65|) have been successful in 
predicting electrostatic properties of biomolecules (I59|) - (|65f) . However, for the complex RNA tertiary struc- 
tures, the line-charge structural model used in the CC theory could be oversimplified (|59r) . Furthermore, 
RNAs with high charge density generally induce high ion concentration and consequently could possibly 
cause ion-ion correlations for multivalent ions in the vicinity of RNA surface. The PB theory ignores such 
potentially important effects for multivalent ion solutions (|66l -l68l). In order to take into account the effects 
of ion correlations and ion-binding ensemble, the tightly bound ion (TBI) model was developed (|66l . 1671): 



see the Supplementary Material for an introduction of the TBI model. Experimental comparisons sug- 
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gest that the TBI model may offer improved predictions on the ion effect of thermodynamic stability of 



DNA/RNA helices/hairpins (|3lM33|) . DNA helix assembly (|69M7l|) . and ion-binding properties and tertiary 
folding stability of RNAs in the presence of Mg 2+ (|72l-l74l). 



In this work, we will employ the newly refined TBI model (|72L 17 3f) to investigate a paradigm system 
with two nucleic acid helices tethered by a loop. We will focus on the effect of conformational constraint 
on the role of ions in RNA folding. The conformational constraint can arise from the macromolecular 
crowding effect. Moreover, the study goes far beyond the previous studies which only involves planar con- 
figurations Specifically, we will consider the conformational fluctuation of helices at the all-atom level 
with 3-dimensional rotations and compare our predictions with the recent experimental data as well as the 
predictions from the PB theory (|44l) . 

Methods 
Structural model 

The model system consists of two helices tethered by a flexible loop; see Fig. 1. In the model, we allow 
for three-dimensional rotation of the helix axes. We note that in a previous study, only symmetric co-planar 
model rotations were allowed (|70|) . As shown in Fig. 1, a configuration of the system can be constructed 
through the following operations. First, the two parallel helices are separated by a axis-axis distance x; 
Second, the two helices are symmetrically rotated around the respective ends (O and O') of the helix axis 
in the axis-axis plane with an angle 6; Third, the two helices are rotated so that the projections of two helix 
axes in the plane perpendicular to line of O-O' are separated by an angle y. The configuration of the system 
can be described by three parameters: the distance x, the angles 6 and y. Here, the parameter y is used to 



produce non-planar configurations. The (DNA) helices here are assumed to adopt the B-form structure (|75|) . 
We use the intervals of (3A, 20°, 20°) to sample the conformational space of (x £ [5A, 56A], 6 e [0, 180°], 
y e [0, 180°]). To reduce the computational time, we randomly select ~600 configurations from the full 3D 
conformation ensemble for calculations. Tests with the different sampling show that our results are quite 
robust. In the structural model, the three parameters (x, 6, y) are used to produce the major portion of 
the full conformational ensemble, including the non-planar/planar and compact/extended conformational 
ensembles. Although this approximation may be valid as we focus on the low -resolution energy landscape 
of the system, it would be useful in further studies to examine the effects of other degrees of freedoms such 
as the spin of helices and asymmetrical rotations of helix axes. 

Free energy 

For a given configuration of the system described by (x,6,y), the total free energy G(x, 8, y) can be 
calculated as 

G(x, 9, y) = G E (x, 6, y) + G loop (x, L) + G cs + G y , ( 1 ) 

where G E (x, 6, y) is the electrostatic free energy for the two helices in an ionic solution. G\ oov (x, L) is the 
free energy of the loop. G cs is the coaxial stacking free energy between the two helices, and Gy is the 



tertiary contact energy as the two helices are in close contact. 

In the study, we calculate G E (x,9,y) using the all-atom TBI model (1721 1731) : see the Supplemental 
Information for a detailed description of the model. We calculate G\ oop (x, L) from 



Gioop — —kBTh\P\ oov {x, L), 



(2) 



where P\ oov {x, L) is the probability of a loop of length L with end-to-end distance x. P\ oov (x, L) can be 



estimated from an approximate analytical expression ( 17 6() : P\ oov (x, L) = ^^jBm ex P ( 4( i -(x/l) 2 ) ) ' wnere ^ 



is the normalization constant (|76M78|) and l p is the persistence length of the loop. In our calculations, in 
order to make direct comparison with the experiment, we assume that the loop is a PEG chain Q and the 



persistence length is 3. 8 A (|44l l70|) . 

The coaxial stacking free energy term G cs is an important energetic factor for nucleic acid helix align- 



ment (28 



791 ). We model this term based on the thermodynamic parameters (|23l l79|) and the RNA structure 



analysis (|28j 



G cs = Go for x < 10A and 9 > 150°; 
= else, 

where Go = Gxx - 6g. Here, is the coaxial stacked base pairs, and Gxx 



(3) 
(4) 

-1.5 kcal/mol for ^ (at 

25°C) in the experimental system (|44|) . Sg(- -1 kcal/mol) is an offset between the coaxial stacking energy 



and the corresponding canonical base pairs (|79|) 



G 3 o is added to account for the possible tertiary contacts (45) 



G 3 o = g y for R cc < 28A; 
= else, 



(5) 
(6) 



where g 3 « is a constant negative energy. In the calculation, we use -\0k B T, -12fc B T, and -\4k B T to investigate 
the effect of the tertiary contact, respectively. We use the condition R cc < 28A, a typical inter-axial distance 
for ion-induced DNA aggregates (|80l— l83h. to simulate a close contact where tertiary contact can occur. Our 
control tests show that the small changes in the criteria around 28A would not cause significant changes in 
our quantitative predictions or qualitative conclusions. 

For the purpose of visualizing the free energy landscape, we will also use the inter-axis angle to 
represent the two angles 9 and y (shown in Fig. 1) 



e 



2 9 2 y \ 
arccos(cos (-)cos (^Oj 



(7) 



to show the electrostatic free energy landscape in 2D (x, 0) plane instead of in 3D (x, 9, y) space. 



Spatial confinement 



As shown in Fig. 1, we use R cc , the distance between the centers of the two axes of the helices, to 
describe the compactness of the system. Furthermore, we use the distance between the outer axial ends 



A 



i?pp' < Rmax, to quantify the spatial constraint of the system. A smaller R max corresponds to a stronger 
spatial confinement due to, for example, macromolecular crowders. 

To quantify the compactness of the system for an ensemble of conformations, we use a Boltzmann- 
weighted averaged value R cc for the distance R cc between the centers of the helix pair (|70l) 

Results and Discussions 



To investigate how the different ionic conditions affect the folding properties of the model system, we 
apply the TBI model to calculate the electrostatic free energy landscape for a wide range of Na + and Mg 2+ 
conditions: [Mg 2+ ] 6 [0.0 ImM, 0.1M] and [Na + ] e [0M, 2M]. To directly simulate the experiments |44|), we 
assume that the system is always immersed in a 16mM Na + background (from the Na-MOPS buffer). For 
convenience, we define the following two specific states for the model system: (i) the random relaxation 
state where the two helices can fluctuate randomly, corresponding to the state with fully neutralized helices 
HQ; (ii) the (electrostatically) folded (collapsed) state with J? cc < 30A flHBS- 

In the following, we will first present the results on the electrostatic free energy landscape. We will then 
discuss the role of ions in the structural collapse (without the possible tertiary contact). Finally, we will 
discuss the effects of spatial confinement and tertiary contacts on the ion role in structural folding. 



Electrostatic free energy landscape 



In Na + solutions. 

Figs. 2A-D show the electrostatic free energy landscape Ge for the different added [Na + ]. At low added 
[Na + ], the helices tend to avoid each other with the largest x and (see Fig. 2A) due to the weak charge 
neutralization. An increased [Na + ] would enhance the ion binding and charge neutralization due to a lower 
entropic penalty for ion-binding, thus causing a weakened inter-helix repulsion. As shown in Fig. 2B, the 
helices fluctuate around less extended conformations with smaller x and values. When [Na + ] becomes 
very high (> 0.3M), the helices become nearly fully neutralized, thus the helices can fluctuate randomly in 
the full conformational space (see Figs. 2CD), except for very compact configurations (very small x and 0) 
where helix-helix and ion-helix exclusions play a dominate role. 



In Mg 2+ solutions. 

Figs. 2E-H show the electrostatic free energy landscapes Ge for the helices with the different added 
[Mg 2+ ]. At low [Mg 2+ ], the helices form extended conformations with large x and 0. With the increase of 
[Mg 2+ ], the two helices switch to the conformations with smaller and x (Figs. 2FG). For high [Mg 2+ ] 
(> 3 mM), the helices tend to adopt the conformations with very small (0 < 40°) and intermediate 
x (32A < x < 47 A); see Fig. 2H. The above transition with increasing [Mg 2+ ] can be attributed to the 
divalent ion-mediated helix-helix interactions (|43l l69l l70|) . When [Mg 2+ ] is very low, background Na + of 
low concentration in buffer only gives weak neutralization to helices, thus the helices repel each other and 



the conformations with largest x and would be favorable. With the increase of [Mg 2+ ], Mg 2+ will com- 
pete with b ackg round Na + , and the self-organization of Mg 2+ could induce (slight) attractive force between 



helices (|69l. r70h. Such attraction would become stronger for the helix-helix near-parallel configurations at 
intermediate separation and at higher [M g 2+ ] (|70r) . As a result, the helices would transform to the confor- 



mations with small and intermediate x (|69|) 



Compared with the results for 2D (co-planar) configurations (Fig. 2 in Ref (1701)). more extended con- 
formations (large 0) are shown in the current 3D free energy landscapes, causing an increase in the con- 
formational entropy of the two helices and weakening the effect of the possible Mg 2+ -mediated helix-helix 
attractive force in the structural collapse. 

Ion-mediated structural collapse: without spatial confinement 

First, we discuss the role of ions in the structural collapse without spatial confinement and tertiary con- 
tact by setting G 3 ° = in Eq. [TJ In the following, we present general features for the system and the 



comparisons with the experimental data (44) 



General features 

As shown in Fig. 3 A, with the increase of [Na + ] to 2M, R cc decreases from a large value to that of the 
random relaxation (fully neutralized) state. Such a trend of R cc is a results of the electrostatic free energy 
landscapes shown in Fig. 2A-D and can be attributed to the aforementioned enhanced Na + -binding at higher 
[Na + ]. Also shown in Fig. 3A is the loop length dependence of the structural compaction at low and high 
[Na + ]. A PEG chain is quite flexible (l p - 3.8A). The loop (length L » l p ) tends to minimize the loop 
entropic loss upon structure formation. Therefore, the loop drives the system to fluctuate around a relatively 
smaller end-to-end distance x since a larger x corresponds to a smaller loop entropy. For example, for PEG 
loops with L = 56 A and 32A, x - 20 A and 15 A, are the most favorable according to Eq. |2J respectively. 
The loop effect is slightly stronger for low [Na + ] where the favorable configurations are extended and the 
effect of the (short) tether can become more important. 

In a solution with added [Mg 2+ ], at low [Mg 2+ ], the favorable conformations are those with large x 
and 0, thus a loop length can contribute to the compaction of the system. In contrast, at a high [Mg 2+ ], 
Mg 2+ induces (slightly) attractive force between helices at small inter-axis angles (see Fig. 2H). When 
the helices are driven by the loop to a smaller x (to gain loop conformation entropy), the low free energy 
conformations (in the small x regime in Figs. 2GH) correspond to a range of 0. As a result, a small x 
does not necessarily lead to compaction of the system (characterized by a small R cc ). We note that a high 
[Mg 2+ ] can induce a state that is slightly more compact than the random relaxation (fully neutralized) state 
yet much looser than the upper limit of the electrostatically folded state as shown in Fig. 3B. 



Comparisons with experiments 

To make direct comparison with the experimental system (|44l) . we assume the loop length to be L = 32A. 
We calculated the Boltzmann-weighted averaged SAXS profiles I(Q) over conformation ensemble for each 
ionic condition. Here Q is the scattering vector which is equal to Ansin{d)l X, where ■& is the Bragg angle 



and A is the wave length of the radiation. 

Figs. 3C&D show the SAXS profiles for the different added [Na + ] and [Mg 2+ ], respectively. As the 
added [Na + ] is increased from to 2M or Mg 2+ is added from to 0.1M, the SAXS profile changes from 
that of the extended state to that of the nearly fully neutralized state. In addition, the profiles indicate 
that the system at high [Mg 2+ ] (~ 0.1M) is slightly more compact than that at high [Na + ] (~ 2M). Since 



the experimental SAXS profiles are in arbitrary unit and with vertical shifts (|44l) . we cannot make direct 
comparisons for the SAXS profiles between the experiments and our predictions. 

Based on the calculated SAXS profiles and compactness R cc , we can estimate the fraction "relaxed" ^ 



for the system by (|44f) 



a (Na + orMg 2+ ) = (1 Q (no added salt) + x O. (high salt), (9) 

where Q. stands for a structural parameter, such as I(Q) or R cc . Q (no added salt) is that in 16mM Na + 
background from the Na-MOPS buffer without added salt, and Q(high salt) is taken as that of the fully- 
neutralized state in our estimation. Based on the equation, we estimate ;f from the calculated I(Q) and R cc 
respectively, and compare our predictions with experiments for the cases with added Na + and Mg 2+ . 

Figs. 3E&F show the fraction "relaxed"^ as a function of added [Na + ] and [Mg 2+ ], respectively. When 
[Na + ] is added, x increases from zero to ~ 0.9 at 2M [Na + ], suggesting that the system is close to the 
random relaxation (fully neutralized) state. The lines estimated from I(Q) and R cc are almost identical, and 
the predictions agree well with the experimental data, as shown in Fig. 3E. 

As shown in Fig. 3F, when Mg 2+ is added, x increases from zero to ~ 0.9 at ~8mM [Mg 2+ ], and may 
slightly exceeds 1 when [Mg 2+ ] > lOmM, suggesting slightly more compact state induced by high [Mg 2+ ]. 
In addition, the predicted x from R cc is slightly larger than that from I(Q) at high [Mg 2+ ] (shown in Fig. 
3F), which is in accordance with the above described relations of R cc and I(Q) with [Na + ] and [Mg 2+ ]. Fig. 
3F also shows that the predicted ^'s from the TBI model are in good accordance with the experimental data, 
while PB theory predicts > 10-fold higher midpoint of [Mg 2+ ] for the Mg 2+ -mediated collapse. This result 
may be attributed to the ignored ion-ion correlations in the PB. The inclusion of ion correlations would 
allow ions to self-organize to form low-energy state, and thus the TBI model predicts a more efficient role 



of Mg 2+ than the PB d67L|72j,|73|). 

It is important to note that although the system has the similar global compactness at high [Na + ] and 
[Mg 2+ ], the free energy landscapes are very different as shown in Fig. 2. 

The PEG loop used in the above calculation is electrically neutral, while a realistic RNA loop is a 
polyanionic chain. Recent experiments suggest that the hinge stiffness can become an energetic barrier for 



RNA folding (I39|) . As shown previously (|70|) . a nucleotide loop may enhance the sharpness of ion-mediated 
structural collapse due to the ion-dependence of the persistence length of the loop (e.g., (0)). 

Ion-mediated folding: with spatial confinement 

In this section, we investigate the influence of spatial restriction as well as a tertiary contact on the ion- 
mediated collapse of the above two-helix system. 
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Electrostatic free energy landscape with spatial confinement. 

With the conformational confinement defined by the limit R max for the distance R PP > between the outer 
ends of the two helices (see Fig. 1), the conformational space of the system is confined. The electrostatic 
free energy landscape of the confined system G £ (7? cc ,7? max ) is given by 

Rp p> <i? max 

G E (R cc ,R max ) = -k B T In J] e- G ^ )/k * T , (10) 

flee 

where R cc is the distance between the two centers of the two helices (see Fig. 1). 

As shown in Fig. 4, the spatial restriction influences the free energy landscape significantly for both Na + 
and Mg 2+ solutions. For large R max , the free energy landscape is close to that without spatial restriction. At 
low [Na + ], the favorable conformations are clustered in the region of large R cc . At high [Na + ] (e.g., 2M), 
the helices can fluctuate nearly in the whole region with all the different R cc values. With the decrease of 
^max, the conformations of the two helices are confined by the restriction and the extended conformations 
are severely limited. This causes the free energy G £ (7? cc ,i? max ) at high [Na + ] to become close to that at low 
[Na + ]. In the limit of strong spatial confinement such as R max < 40A, the free energy G £ (/? cc ,7? max ) at the 
different Na + concentrations can become nearly identical; see Figs. 4A-D in the small R max regime. 

For Mg 2+ solutions, the spatial confinement effect on G£(i? cc ,7? max ) is similar to that for Na + solutions, 
except for the much more efficient role of the Mg 2+ ions in promoting the collapse of the system. The 
electrostatic free energy landscapes show that with the spatial confinement, the extended conformations are 
much more significantly restricted than the compact conformations, causing the destabilization of the ex- 
tended state. In the ion-promoted collapse, the spatial confinement would effectively enhance the ion effect 
in folding. 



Ion-mediated folding: effects of spatial confinement and tertiary contacts 

To investigate the role of tertiary contacts in ion-induced structural collapse, we add a negative constant 
energy to model the effect of the tertiary contacts. Specifically, we choose three values -10k B T, -\2k B T, and 
-14£ B Tforg 3 ° (see Eq.©. 

Helix-helix attraction arising from the tertiary contacts would cause a shift in the conformational dis- 
tribution toward the compact state and lower the ion concentration required to induce the collapse of the 
structure. Furthermore, the tertiary contact results in a much sharper ion-induced structural collapse (than 
the tertiary contact-free case). The effect is more pronounced for stronger tertiary contacts. As shown in 
Figs. 5ACE for Na + , the midpoint for the folding transition occurs at [Na + ]~ 0.5M for g 3 = = -\0k B T 
and [Na + ]~ 0.08M for g^ = -l4k B T, respectively, and the stronger tertiary contacts cause a slightly more 
compact state. As shown in Figs. 5BDF, the midpoints of of Mg 2+ -dependent folding are ~0.8mM and 
~0.2mM for gy = -\0k B T and gi° = -l4k B T, respectively. Again, we note that much lower [Mg 2+ ] (~ 
mM) is required to cause folding than [Na + ] (~ M), which is in accordance with the experimental data 

To further examine the effect of the spatial confinement in the presence of tertiary contacts, we calculated 
R cc with a given i? max constraint (the maximum distance of P-P' shown in Fig. 1). As shown in Figs. 5A-F, 
the spatial confinement significantly enhances the ion efficiency in mediating the structural collapse/folding 



in both Na + and Mg 2+ solutions. For example, for gy = -\0k B T, with the decrease of R max from 120A 
to 70A, the [Na + ] and [Mg 2+ ] at the midpoint of the folding transition decreases from ~0.5M to ~ 0.1M 
and from ~lmM to ~0.3mM, respectively. For a stronger tertiary contact, e.g., gy = -\4k B T, the spatial 
confinement would further enhance the ion efficiency for both of Na + and Mg 2+ , as shown in Figs. 5A-F. 
Such higher efficiency of Mg 2+ than Na + is attributed to the higher valence of Mg 2+ which corresponds to 



the less entropy for ion binding and stronger ion-ion correlation, causing higher Mg 2+ binding affinity (|72l) 
and consequently higher efficiency in stabilizing RNA (|73r) . 

The above spatial confinement-promoted ion (Na + and Mg 2+ ) efficiency in mediating structural folding 
is in accordance with the recent exp eriment which showed that a lower Mg 2+ is required to fold Azoarcus ri- 



bozyme in PEG environments (|19l) . The spatial confinement can be caused by molecular crowding or other 
in vivo effects. The spatial confinement can destabilize the low- salt (extended) state by reducing the con- 
formational space and consequently decrease the electrostatic free energy difference between the compact 



conformations and the (restricted) extended conformations (|52r) . In addition, the tertiary contacts would 
further stabilize the compact state. These two effects combined together cause a reduced ion concentration 
required to induce the folding transition. Therefore, Na + and Mg 2+ are more efficient in RNA folding in a 
crowded environment. 

Conclusions and Discussions 



In this work, we employed the newly refined TBI model (|72L |73|) to investigate how ions cause com- 
paction for a system of loop-tethered helices. The main advantage of the TBI model is its ability to account 
for the possible ion correlation and fluctuation effects. The present study distinguishes from the previous 
studies (|70l) on similar systems in several novel aspects. 

First, the present study is based on the 3D (instead of planar) rotational degrees of freedom in the 
conformational sampling of the system. In contrast to co-planar (2D) helix orientations, the 3D random 
orientation of the helices would reduce the (ensemble averaged) ion correlation effect and the associated 
force between the helices. 

Second, the main focus of the present study is the influence of the spatial confinement and the ter- 
tiary contacts upon the ion effect and the comparisons with the the PB-based predictions and the recent 
experimental results. The study leads to the following main conclusions: 

1. Na + can induce a transition from an ensemble of extended conformations to an ensemble of random 
relaxed conformations with a flat energy landscape. In contrast, Mg 2+ can induce a transition from 
an extended state to a conformation ensemble that is slightly more ordered and slightly more com- 
pact than the random relaxation conformational ensemble. However, such a Mg 2+ -induced state is 
far looser than the (electrostatically) folded state. The predictions from the TBI model are in accor- 
dance with the experimental data for both Na + and Mg 2+ solutions, while the PB theory significantly 
underestimates the Mg 2+ efficiency. 

2. A tertiary contact would make the ion (Na + and Mg 2+ )-dependent folding transition much sharper 
than the tertiary contact-free case. Furthermore, stronger tertiary contacts can fold RNAs at lower ion 



Q 



concentrations for both of Na + and Mg 



3. Spatial confinement (e.g., due to macromolecular crowding) significantly promotes the ion efficiency 
in mediating RNA folding for both Na + and Mg 2+ . The effect of the spatial confinement can be 
attributed to the decreased electrostatic free energy difference between the compact ensemble and the 
(restricted) extended ensemble. 

The current theory gives overall good agreement with the experiment results, suggesting an overall reli- 
able analysis and predictions for the ion-induced compaction and the effects of spatial confinement. Further 
investigation of the problem requires improvements of the model in several aspects. First, in the structural 
modeling, although we have taken in account the non-planar conformations, we have ignored conformations 
generated from the spin of the helices and the nonsymmetric rotations of the helix axes. Including these 
conformations may cause changes in the free energy landscape, though a detailed study for such changes 
requires exceedingly high computational cost. Second, the current form of the TBI model cannot treat the 
sequence preference in ion-binding (f7ll l8~5|) . The specific binding could make important contributions to 
the ion effect, although for the PEG-tethered helix system here, the contribution of the specific binding may 
be weak. Finally, we have neglected the interference from the helix on the loop conformations except for 
the constraint of the end-end distance x. The effect of the loop-helix interference may not be significant for 
a PEG loop. However, such a effect can be important for a polynucleotide loop which is highly charged. 
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FIGURE CAPTION 



FIGURE 1 The conformation of two (DNA) helices tethered by a loop can be characterized by three struc- 
tural parameters: the end-to-end distance x (O-O'), the angles 6 and y illustrated in the figure. The 
freedom of y (rotation of the axis into the paper) is used to produce non-planar configurations. The 
P-P' distance < R max represents the spatial restriction for the two helices due to, e.g., the macro- 
molecular crowding effect. R cc between the centers of the helices is used to quantify the compactness 
of the model. Tertiary contact energy can be added for conformations with a small R cc . The DNA 
helices are constructed with the software X3DNA (jvHl) and are displayed with the software PyMol 
(http://py mol.sourceforge.net/| ). 

FIGURE 2 The electrostatic free energy landscapes G E (x, 6,y) as a function of [Na + ] (A-D) and [Mg 2+ ] 
(E-F). Here x is the end-to-end distance and is the angle between the two helix axes. is deter- 
mined by the (0,y) angles (see Eq. [7]). The color scale represents the energy difference between the 
configuration (x, 0) and the state with the minimum free energy. From to 5 k B T , the color changes 
ly from Red to Blue. Note that the system is in a 16mM Na + background due to the Na-MOPS 



gradua 
buffer d44h. 



FIGURE 3 Ion-mediated structural collapse in the absence of the tertiary contact and spatial confinement. 
(A, B) Ion concentration-dependence of the compactness R cc (the mean distance between the he- 
lix centers) as a function of Na + (A) and Mg 2+ (B). Dashed lines denote the random relaxation 
(fully-neutralized) state and the dotted lines denote the upper limit of the electrostatically folded 



state (R cc ~ 30A) (j69|, LZO, l80j4S2b- (C,D) The SAXS profiles calculated for the different added Na + 
(C) and Mg 2+ (D). Lines are the reference SAXS profiles. Dashed lines: the upper limit of the elec- 
trostatically folded state (R cc - 30A); Solid lines: the random relaxation (fully-neutralized) state; 
Dotted lines: the extended state (with no added salts). The SAXS profiles are calculated from the 



software CRYSOL (84) with the statistical weight determined from the free energy G(x, 6,y). (E,F) 



The "relaxed" fraction^ as a function of [Na + ] (E) and [Mg 2+ ] (F), respectively; see Eq. |9] Symbols: 



the experimental data by Bai et al (|44l) : Solid lines: the calculated results with the SAXS profiles I(Q) 
determined from the TBI model; Dotted lines: calculated with the compactness R cc determined from 
the TBI model; Dashed lines: calculated with the SAXS profiles determined from the PB theory by 



Bai et al (|44l) . The "relaxed" state is taken as the fu 
16mM Na + background due to the Na-MOPS buffer (|44|) 



ly neutralized state. Note that the system is in 



FIGURE 4 The electrostatic free energy landscapes G E (R cc ,R max ), where R max represents the spatial re- 
striction for the helices (e.g., due to the macromolecular crowding); see Eq. \\0\ The color scale 
denotes the energy difference between G E (R cc ,R max ) and the minimum value for each R max : AG E = 
G E (R cc ,R max ) - G E (R cc ,R max ) m i n . From to 5k B T, the color changes uniformly from Red to Blue. 

FIGURE 5 Ion-mediated structural folding in the presence of tertiary contact and spatial restriction. (A,B) 
The compactness R cc of the loop-tethered helices as functions of [Na + ] (A) and [Mg 2+ ] (B) for the 
different R max values (from the right to the left: 120A, 82A, 78A, and 70A) and the different tertiary 



contact energy g 3 » (solid lines: -lOk B T, and dashed lines: -14 k B T). (C,D) The folded fraction /folded as 
a function of [Na + ] (C) and [Mg 2+ ] (D) for the different R max values (from the right to the left: 120A, 
82A, 78A and 70A) and the different tertiary contact energy g 3 ° (solid lines: -10k B T, and dashed 
lines: -14 k B T). (E,F) The ion concentration midpoints for the folding transition as a function of R m3X 
and the different tertiary contact energies gy. From the top to the bottom: g y =-\0k B T , -12k B T, and 
-\Ak B T, respectively. 
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